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Dark Discrete Gauge Symmetries 
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We investigate scenarios in which dark matter is stabilized by an abelian Zn discrete gauge 
symmetry. Models are surveyed according to symmetries and matter content. Multi-component dark 
matter arises when N is not prime and Zn contains one or more subgroups. The dark sector interacts 
with the visible sector through the renormalizable kinetic mixing and Higgs portal operators, and 
we highlight the basic phenomenology in these scenarios. In particular, multiple species of dark 
matter can lead to an unconventional nuclear recoil spectrum in direct detection experiments, while 
the presence of new light states in the dark sector can dramatically affect the decays of the Higgs 
at the Tevatron and LHC, thus providing a window into the gauge origin of the stability of dark 
matter. 
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There is compelling empirical evidence for the exis- 
tence of dark matter (DM) in the Universe [Q . A simple 
and attractive possibility for DM is a new elementary 
particle that is stable on cosmological times scales. This 
stability strongly suggests the existence of a new sym- 
metry in the dark sector, and understanding the precise 
nature of this symmetry is of fundamental importance in 
our quest for a theory of DM. 

While there are many possibilities for symmetries that 
stabilize DM, by far the most common example invoked 
in the literature is a discrete Z2 symmetry, e.g. R-parity 
in Minimal Supersymmetric Standard Model (MSSM). 
Generic parity symmetries are clearly a useful tool in 
building models of DM, but they often appear ad hoc. 
Indeed, discrete symmetries are not an organizing princi- 
ple of the Standard Model (SM), and in general quantum 
field theories suffer no theoretical inconsistency once dis- 
crete symmetries are abandoned. It has also been argued 
that global discrete symmetries will be violated by quan- 
tum gravitational effects Q. It is therefore reasonable 
to ask if discrete symmetries may find their origin as a 
consequence of a more fundamental principle that may 
ultimately be responsible for the stability of DM. 

Gauge symmetry is just such a principle that may pro- 
vide the eventual rationale for discrete symmetries. A 
so-called discrete gauge symmetry emerges as a remnant 
of a spontaneously broken gauge symmetry g]. These 
symmetries are exactly conserved owing to their gauge 
origin, and thus the lightest state having a nontrivial 
charge under the discrete symmetry is absolutely stable. 
In the past, discrete gauge symmetries have been used 
in a variety of contexts, but most notably to explain or 
find an alternative to i?-parity in the MSSM ||. Very few 
studies have been solely motivated by explaining DM sta- 
bility, but see Ref. [|j for a recent example based on a 
Zi remnant from a nonabelian SU(2) gauge symmetry. 
See Ref. B for earlier work on dark matter models with 
nonabelian discrete gauge symmetries. 

In this paper we will investigate models of hidden 
sector DM based on a spontaneously broken U(l)rj 
gauge symmetry which preserves a discrete Zn subgroup. 
While U(1)_d is the minimal choice for the dark sector 



gauge group, it also affords the possibility of renormal- 
izable interactions of DM with the SM via the kinetic 
mixing M and Higgs portal operators. New phenom- 
ena may occur for large Zn symmetries. This includes 
multiple stable DM candidates, new direct detection sig- 
natures, and novel decays of the SM Higgs boson leading 
to final states with multiple leptons, jets, and missing 
energy. These signatures provide a means to probe the 
gauge origin of DM stability. 

Discrete Gauge Symmetry. Let us begin by review- 
ing the prototype setup for a gauged Zn discrete symme- 
try [[2] . Consider a U(1)d gauge theory containing a dark 
Higgs field 4>n with integer charge N and scalar matter 
field xi with charge 1, where the notation is that the 
subscript indicates the charge of the field under U(1)d- 
In addition to the usual self-conjugate terms the gauge 
symmetry allows the following term in the Lagrangian: 



AC ex <P nX i + h.c. 



(1) 



We assume the scalar potential is constructed so that 
the dark Higgs gets a vacuum expectation value (vev). 
This spontaneously breaks the gauge symmetry, but the 
Lagrangian preserves a discrete Zn subgroup of U(1)d 
under which the matter fields transform as 
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e 2 ™ k / N X i, 



k = 0,l,2,...N- I. 



(2) 



In general, as long as there are no additional sources of 
gauge symmetry breaking and all gauge anomalies are 
canceled with appropriate matter content, the remnant 
discrete symmetry will be exactly conserved. The case 
N = 2 corresponds to the usual parity so often encoun- 
tered in models of DM, but in this construction there is 
nothing particularly special about the choice of Z2 , as it is 
simply a consequence of the relative charges of the Higgs 
and matter fields in the dark sector which are a priori 
undetermined. Note that depending on the field content 
and associated gauge and Lorentz charges, there may be 
no renormalizable terms analogous to Eq. (pi) present. 
In this case there will be accidental global symmetries 
which also may stabilize DM, at least at the renormal- 
izable level. Of course these accidental symmetries may 



be violated at the nonrenormalizable level, in contrast to 
the discrete gauge symmetry. 

Models. We now give a survey of DM models based on 
an abelian gauged Zn symmetry, generalizing Eq. (Q) to 
include fcrmions, multiple matter fields and chiral mat- 
ter. Models are characterized by the field content and 
corresponding charges. The full Lagrangian contains ki- 
netic terms for the U(1)d gauge boson V^ and the dark 
Higgs boson </>jv, as well as an appropriate scalar poten- 
tial such that 4>n acquires a vev (4>n) = v 1 /V2. This vev 
breaks the gauge symmetry leaving a massive vector and 
a physical dark Higgs fluctuation h! . In addition there are 
gauge kinetic mixing and Higgs portal terms that allow 
interactions with the SM; we will discuss these in detail 
below. Finally, there will be model dependent terms re- 
lated to the additional matter fields which comprise the 
DM sector. We begin by examining models with minimal 
matter content, i.e. one dark Higgs and one matter field, 
leading to three possibilities at the renormalizable level: 
1) a Z2 scalar, 2) a Z2 Dirac fermion, and 3) a Z3 scalar. 
We then turn to multi-field and chiral models. 

Zi scalar: The Zi scalar DM model contains a charge 
2 dark Higgs field 4>2 and a charge 1 scalar matter field 
Xi- The gauge symmetry permits the term (|l|) to appear 
in the Lagrangian: 
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4xiXi + h.c. 
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where a possible phase in the coupling A may be re- 
moved by rotating \i- A novel aspect of Z2 mod- 
els is that Eq. (||) induces a mass splitting between 
the real and imaginary components of the matter field, 
Xi = (S + iP)/y2, after U(1)d is spontaneously broken. 
Thus, the physical masses of the S and P scalars are 
m 2 s P = rriQ =p Am 2 where TOq includes universal contri- 
butions from a bare mass operator X1X1 an d Higgs por- 
tal operators O^Xxlxi), an d {H^H)(x\xi) following 
gauge symmetry breaking, while Am 2 = Aw' comes from 
Eq. (|^). According to Eq. (||) the components of the 
multiplet are Z 2 odd, S —> —S, P —> —P, so that for 
negative A, S will be the stable DM candidate. 

The splitting of the individual components S, P leads 
to an off-diagonal interaction coming from the gauge ki- 
netic term of the matter field: 



goV^P^S, 



(4) 



where qd is the gauge coupling of U(1)d- One well- 
known consequence of this operator is that scattering 
with nuclei mediated by vectors V^ becomes inelastic [|| . 
Another important effect is that this interaction will in- 
evitably allow the heavier Zi partner P to decay to the 
DM plus light SM fermions which will have important 
implications for collider physics, as we will discuss in de- 
tail later. 



Zi fermion. Consider next the minimal theory of 
fermionic DM stabilized by a discrete Z 2 gauge symme- 
try. The model contains a charge 2 dark Higgs (f>2 , and a 
vectorlike pair of Weyl fermions ipi and £_i with charges 
+1 and -1, respectively Besides the kinetic terms and 
Dirac mass term, C D — m_D(^i£_i+h.c.), the gauge sym- 
metry permits the following terms analogous to Eq. (pi): 
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Like the Z2 scalar above, there is a splitting between the 
physical mass eigenstates states, which can be written 
in terms of Majorana fcrmions \IJ±, with masses m± ~ 
m£)±(rn.L-|-77i.R)/2. Here we have defined n%L,R = Xl,rv' 
and made the assumption mo 3> rriL,R for simplicity. 
There is discrete Z2 gauge symmetry that remains under 
which both ^>± are odd, ensuring the stability of the 
lightest particle. In the limit where Xl = Xr, the gauge 
interaction becomes 



ignVf, *+7^*_. 



(G) 



We see that, at least structurally, the Zi fermion has 
many similar features to the Z2 scalar discussed above, 
and in particular the gauge interaction is off-diagonal. 
We note that both the Z^ scalar and Z2 fermion mod- 
els were considered in the past as models of inelastic 
DM ||, g] , and a similar model with two Dirac fermions 
was examined in Rcf. pi. 

Z% scalar. The final model with minimal field content 
is the the Z% scalar DM model. This model contains a 
charge 3 Higgs field ^3 and a charge 1 scalar matter field 
Xi- The gauge symmetries allow the additional term (Fy) 



\/2 t 
A£ = — -y-X ^3X1X1X1 + h.c. 



(7) 



Unlike the Z2 models considered above, there is no mass 
splitting between the real and imaginary components, 
which clearly must be the case since Z 3 charged fields 
are necessarily complex. After symmetry breaking, the 
interactions from (M) become 

~ 3! Xv ' XlXlXl ~ ^ Xh ' x ^XiXi + h.c. (8) 

The interaction with the dark Higgs h' leads to the possi- 
bility of 'semi-annihilation', X1X1 ~~ ^ Xi^-'j an( i may have 
interesting consequences for the thermal history of xi, as 
emphasized recently in |Q|. Note that the interactions 
in (||) being cubic in xi, do not allow for a tree- level DM- 
nucleus scattering. Direct detection is of course still pos- 
sible through other portal interactions. For other models 
of Z 3 DM, see Rcf. || 

Multi-field Zn models. With a single matter field, as 
considered in the models above, the renormalizable pos- 
sibilities for discrete gauge symmetries are limited to Z2 



and Z3. Once we allow for more than one matter field in 
the dark sector, renormalizable Zn models are possible. 
Here we will consider several interesting features of these 
models, illustrating them with examples. 

Generically, if N is a prime number then barring acci- 
dental global symmetries there will be only one DM can- 
didate: the lightest particle with nontrivial Zjq charge 1 . 
However, if N is not prime, but rather composite (has 
divisors besides 1 and N), then the Lagrangian will pre- 
serve additional discrete symmetries that are subgroups 
of Zn, allowing for the possibility of multiple stable 
species and thus multiple DM candidates |il| . Note that 
multi-component DM has been considered for a variety of 
reasons in the past; see Ref. H]. To illustrate, consider 
the simplest case in which this occurs, namely a Z4 sym- 
metry, which has Z 2 as a subgroup. A multi-component 
DM model may be constructed with the following field 
content: a charge 4 dark Higgs </>4 and scalar matter fields 
Xi and xi with charges 2 and 1. The Lagrangian contains 
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-Ai</>4X2X2 - A 2 04X2XiXi - A3X2X1X1 + h-' 



(9) 

After symmetry breaking, the first term in Eq. (g) leads 
to a mass splitting for the real and imaginary components 
S,P of the field X2, analogous to the Z 2 scalar model 
considered previously in Eq. (g), and we take Ai negative 
so ms < rap. According to Eq. (g) the theory has a 
descendant Z4 discrete symmetry under which both \i 
and Xi are charged, as well as a Z2 symmetry under 
which S, P are even and xi is odd. The transformation 
properties are summarized in Table |. It is clear that Xi > 
being the only particle with a nontrivial Z2 charge, is 
stable and a DM candidate. However, if S is lighter than 
twice the mass of xi it will also be stable by virtue of the 
Z4 symmetry. Hence, in this case there are potentially 
two species of DM, though what contribution each makes 
to the cosmological DM depends on other considerations 
that we will examine later. 

As a slightly more complicated example, consider a Zq 
model with a charge 6 dark Higgs 4>§ and three scalar 
matter fields xii X2, and X3 with charges 1, 2, and 3. 
The Lagrangian contains 

AC = -Ai4x3X3 - A 2 </4x3X2Xi ~ A 3 ^X2X2X2 

-A4X3X2X1 - A5X3X1X1X1 - A6X2X1X1 + h.c. 

(10) 

In this case, we can decompose Zq = Z3 x Z 2 , so that 
there are two distinct cyclic symmetries. The transfor- 
mation properties of the matter fields under these sym- 
metries are displayed in Table El We see again that there 



1 Additional heavier particles with Zjv charge can be stable if all 
possible decay modes are forbidden by kinematics. 
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TABLE I: Transformation properties of matter fields in the 
Z4 (left) and Ze, (right) multi-field models. Charges qi for a 
field \i are defined through their transformation under Zm'- 
\i ->■ exp(^ ■ qi) X z- 



are potentially two species of DM. For example, if xi is 
the heaviest matter field, then X2 and X3 will be stable 
due to their nontrivial charges under Z3 and Z 2 , respec- 
tively. Note that X3 will be split into its real and imag- 
inary components due to dark symmetry breaking and 
the lightest of these will be stable. 



From the first example based on the Z4 model, we learn 
that a Z p m symmetry with p a prime number and m a 
natural number allows for potentially m stable states, 
while the example with the Zq symmetry shows us that 
depending on N, Zn may be decomposed into a direct 
product of smaller groups [13l . We may surmise that for 
a general N = p^^p™ 2 ■ ■ ■ P™ fc > where pi is prime (pi =/= pj 
for i =/= j) and rrii is natural, Zn is decomposed into the 
product group 
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Pi 



x Z, 



l>2 



2 X 



(11) 



There will be at most mi + vn 2 + • • • + mk stable species 
for a given Zn symmetry, though the actual number de- 
pends on the field content and the spectrum. A simple 
recipe to obtain the maximum number of stable states 



is as follows: for each prime pi in (11), add fields with 
U(l)i5 charges Q = {N/p^N/pf, ■ ■ ■ N/p™'} and order 
the spectrum m N / Pt < pi m N / p 2 < ■ ■ ■ < p™ H ~ m N ^ p ^i . 
Besides multi-component DM, another important ef- 
fect which can occur in multi-field Zm models is mass 
mixing induced by dark symmetry breaking. To illus- 
trate, consider a model with a charge N dark Higgs, and 
scalar matter fields XQ an d XQ+n with charges Q and 
Q + N. Gauge symmetry allows the following term: 



A£ = 



~V2 



4>nXqXQ+n +h.c, 



(12) 



which induces a mass mixing between the matter fields 
Xq and XQ+N following dark symmetry breaking. This 
can be diagonalized with an orthogonal rotation leading 
to mass eigenstates Xa,b- Consider for simplicity the limit 
in which the bare masses are equal, ttiq = ttiq +n = m 2 , 



and a small mass splitting induced by the operator in ( 112 ) 
given by Am 2 = Xv <C m. The mixing angle in this case 
is 45°, and the masses are split by a small amount S ~ 
Am 2 /M . Besides a diagonal gauge interaction between 
the charged matter fields, this rotation induces an off- 
diagonal gauge interaction: 



■ t/ ( t a r y + iV 



Q + 7V/2 iV/2 

q + n/2 



X.a 

Xb 



(13) 



These interactions can lead to important consequences 
for direct detection experiments, as DM may undergo 
both elastic and inelastic scattering with nuclei, leading 
to a distinctive recoil spectrum. Moreover, off-diagonal 
terms allow for novel decay modes of Zn partners at 
colliders. We will explore these issues below. 

One final common occurrence in Zn models is the ap- 
pearance of nontrivial interactions between matter par- 
ticles. Cubic and quartic vertices are allowed by gauge 
invariance so long as the corresponding charges cancel, 
and these can lead to important effects in DM physics. 
For example, one may have a Yukawa-type interaction 
between scalars \i an d fermion matter fields ipi, 



AC D -Xxitpjipk + h.c. 



(14) 



so long as the charges satisfy i + j + k = 0. This interac- 
tion, along with similar cubic and quartic scalar matter 
interactions, can manifest in DM annihilation, affecting 
the relic abundance and indirect detection prospects, and 
lead to cascade decays of Zn partners given a production 
mechanism at collider experiments. 

Chiral models. Lastly, we will consider models in which 
matter consists of a chiral set of fermions. This means the 
fermion content is chosen so that 1) there are no gauge 
anomalies, and 2) a vectorlike mass term is forbidden 
by gauge symmetry. Indeed, it is a striking fact that 
each generation of SM matter forms a chiral set given 
the gauge group SU(3) c xSU(2)ixU(l)y, and so it seems 
quite plausible that DM could also be chiral. 

Given that the dark fermions are charged only under 
U(1)d, there are two anomaly conditions to satisfy: the 
U(1)5j anomaly and the U(l)zj-gravitational anomaly. 
To forbid mass terms, any two charges must not sum 
to zero. For fermions with charges Qi, these conditions 
are simply 



T.QI- 


= 


Y,Qi- 


= 



Qj 7^ 0, for all i,j 



(15) 



Rcf . [n5l analyzed quite generally the problem of finding 
fermion charges satisfying the conditions (ITq) and pre- 
sented many examples of chiral sets. While there are 
clearly a vast number of possibilities for models of chiral 



DM, we will be content here to illustrate a few common 
features in a simple model. Consider the following chiral 
set: 



^1,-03,^4,^-6^7 



(16) 



where i — 1,2. This set contains 7 'flavors' of Weyl 
fermions, and it is straightforward to check that the con- 
ditions ( |l5| ) are satisfied. Besides the fermions in the chi- 
ral set, we will need at least one dark Higgs to break the 
gauge symmetry and to generate masses for the fermions. 
We choose a charge 7 Higgs field 07 to accomplish this 
task. This means that ultimately there will be a discrete 
Z-j symmetry that stabilizes the lightest fermion mass 
eigenstate, but as we will see there may also be acciden- 
tal global flavor symmetries which can stabilize multiple 
candidates. 

The Higgs </>7 allows for 03 , 04 to marry through the 
Yukawa operator <fi 7 ip3ip4 + h.c, to form a Dirac fermion. 

Furthermore, - 0!_ 1 ,V'-6 a ^ so are marr ied by the Higgs 
through the following Yukawa interaction: 
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-A^r^V-ie + h.c. 



(17) 



where X 3 is a general matrix. This Yukawa mixing is 
reminiscent of what occurs in the SM, but we should 
remember that these fermions are not part of separate 
'generations'. The Yukawa matrix can be diagonalized 
by separate unitary transformations on ip I _ l and ip J _ 6 and 
these unitary matrices are not observable with the field 
content we have so far specified. In particular, there will 
be no dark flavor change induced by the gauge interac- 
tions, as there is an analogue of the GIM mechanism 
operating in this model. 

What about the remaining fermion 07 in the chiral set 
(Eo)? This fermion is quite similar to the neutrinos of the 
SM as no tree-level mass may be generated with the mat- 
ter and Higgs content considered so far. How may we gen- 
erate mass for this fermion? First, let us form the U(1)d 
neutral operator 4yip7 having dimension 5/2 and spinor 
Lorentz index. This operator, being analogous to HL 
of the SM, immediately suggests mechanisms for mass 
generation. A Majorana mass term may arise from the 
dimension five operator (<0^07) 2 , or we may add a singlet 
fermion £0 and form a Dirac mass 0707^0 + h.c. As the 
singlet may also have a bare mass, a see-saw mechanism 
is potentially operative. If U(l)u is broken around the 
weak scale, and the dimension 5 Majorana mass operator 
or the singlet mass is generated at a higher scale, then 
typically we are led to predict a much lighter fermionic 
state in the spectrum. This state is potentially stable and 
may be difficult to deplete in the early universe, having 
few annihilation channels which are unsuppressed. There 
are many ways out of this. For example, the singlet £0 
may couple to the SM neutrinos, breaking the accidental 
flavor symmetry and allowing the state to decay. We may 
add new dark Higgs fields to make the neutrino-like state 



heavy, one example being a charge 14 Higgs which would 
not disturb the subgroup discrete Z7 symmetry. Alter- 
natively, a very light or massless state in the dark sector 
could manifest itself as a new component to the cosmo- 
logical dark radiation and have interesting signatures in 
its own right. 

Finally, let us elaborate on the symmetries of this 
model. By construction, there is a Z7 discrete symme- 
try which will forbid the lightest state with nontrivial 
charge from decaying. However, there are additional ac- 
cidental flavor symmetries present which lead to multiple 
stable states. Let us trace the origin of these symme- 
tries. First, if we turn off all interactions, there is a large 
U(7) flavor symmetry, which is broken to U(2) 2 xU(l) 3 
by the gauge interactions. The Yukawa terms further 
break the flavor symmetry to U(l) 4 leading to poten- 
tially 4 stable mass eigenstates. The U(l) flavor sym- 
metry associated with the neutrino-like state ipj may 
be broken further to a Z2 symmetry by the inclusion 
of a dimension 5 Majorana mass operator, or to a di- 
agonal subgroup of the combined U(l) 2 symmetry with 
the inclusion of the singlet £0 . Finally, higher dimension 
operators will generically further break these accidental 
symmetries. For example, integrating out a charged 2 
scalar field (which does not condense) generates oper- 
ators like (V>_iV'-i)(V'-i'03), (ip4ip-6)(^-i4>-i), etc., al- 
lowing some of the would-be stable states to decay. How- 
ever, the discrete Z7 symmetry is exactly conserved and 
will ensure the stability of the lightest charged state, thus 
providing a DM candidate. 

Phenomenology. We will now describe basic aspects 
of the phenomenology in this class of models. Detailed 
scans of the parameter spaces of specific models is beyond 
the scope of the present work. Rather, we will focus on 
the novel implications of a larger Z^ discrete symme- 
try and the impact of additional states in the dark sec- 
tor such as Zjf matter as well as dark gauge and Higgs 
bosons. The first question we must address is the fol- 
lowing: how can the dark sector communicate with the 
SM? 

Portals. Given that we are considering models in which 
DM is in a hidden sector (DM carries no charge under 
the SM gauge group, and likewise the SM fields carry no 
charge under U(1)d), if interactions are to exist we ex- 
pect that at high scales there exist heavy states charged 
under both sectors. Integrating these out, we generate all 
higher-dimensional operators allowed by the gauge sym- 
metries of the theory. We may also generate 'portal' op- 
erators Jl6[ : renormalizable, gauge- invariant operators 
that connect the SM to the dark sector. In particular, 
a dark sector based on a broken V(1)d gauge symme- 
try affords the opportunity to interact through both the 
kinetic mixing B and Higgs portal M: 



where of course the final term is present only for scalar 
matter fields at the renormalizable level. Note that the 
combination of vector and Higgs portals was first studied 
in Ref. [17]]. Besides kinetic mixing between vectors, sym- 
metry breaking induces mass mixing between the SM and 
dark Higgs bosons. We will limit ourselves to the case of 
small mixing arising from portals, so that it suffices to 
treat the terms in Eq. (nm as interactions. 

Alternatively, one can diagonalize the Lagrangian and 
work directly with the physical states. The main effect 
of the kinetic mixing is that SM matter picks up a small 
charge under U(1)d proportional to the kinetic mixing 
strength k: 

C D kV^-CwJ^m + s w (l - ml/mfr)- 1 J£\. (19) 

The Higgs portal mass mixing may be treated by first ex- 
panding around the vacuum, v — > v + h and v' — > v' + h' , 
followed by diagonalizing the system with an orthogo- 
nal rotation by angle 8^ ~ 2\ ( j ) vv' /m^, valid in the limit 
of m^i -C m/j. Diagonalization induces interactions be- 
tween the SM Higgs and the dark sector, as well as the 
dark Higgs with the SM, 



CD0h(h'J h -hJ h ' 



(20) 
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where Jh,h' are the SM and dark currents coupling to the 
Higgs bosons. Note that for a dark symmetry breaking 
scale well below the weak scale, there can be an issue with 
technical naturalness since the SM Higgs vev will tend to 
raise the dark Higgs vev for large Higgs portal coupling 
A0 in Eq. (|l8|). In some cases we will be interested in 
couplings of order the bottom quark Yukawa, so that 
e.g. a symmetry breaking scale in the dark sector of 
v' ~ 10 GeV requires only a mild order one cancellation 
between different terms in the scalar potential. Scales 
much lighter than this require smaller values of A^ for 
technical naturalness. Similar comments apply to the 
mass of scalar DM when the Higgs portal coupling A x in 
Eq. (§J) is sizable. 

Existing constraints. Since the dark sector carries no 
SM charges, new states can in principle be much lighter 
than the weak scale and still avoid constraints coming 
from direct production and precision measurements. The 
constraints that do exist depend on the strength of the 
portals and the precise decay patterns of dark states. 
We now give an overview of the possible constraints for 
various mass scales in the dark sector. 

Let us first discuss constraints on the kinetic mix- 
ing portal. If the dark gauge boson mass is above 
the weak scale, there are mild constraints on the ki- 
netic mixing parameter n coming from electroweak pre- 
cision tests, k < 10 _1 — 10~ 2 p8[ . Vectors with masses 
(9(1 GeV) < m v < 0(100 GeV) can be constrained fur- 
ther from a variety of e + e~ collider data, and give a 
rough bound of k < few x 10~ 2 |l9f. For sub-GeV vec- 
tors, additional constraints arise from corrections to the 



anomalous magnetic moment of the electron and muon, 
extending k down to < few x 10~ 3 |M. There are also 
model-dependent constraints for light vectors, assuming 
the dominant decay mode of V occurs through the kinetic 
mixing portal back to SM states, arising from searches at 
B and other meson factories as well as fixed target exper- 
iments [B4l-p7| . In particular the fixed target experiments 
are able to exclude a large portion of the parameter space 
for vector masses below a few hundred MeV, with con- 
straints on k ranging from 10~ 4 — 1CP 8 fl26L E7j . There has 
been a great deal of effort invested in the last two years to 
constrain or search for new light gauge bosons decaying 
to leptons, primarily motivated by a possible DM con- 
nection [£l| to cosmic ray anomalies [£2], E3|. There are 
the additional possibilities of detecting or constraining 
light gauge bosons decaying to leptons through searches 
of the existing B-factory data sets [M and new fixed 
target experiments pq| . 

Constraints on the Higgs portal are much weaker and 
only model-dependent bounds are possible. Very large 
O(l) mixing angles Oh can affect electroweak precision 
observables J44{ , but for a sub- weak scale dark Higgs h', 
a stronger constraint comes from direct production at the 
LEP experiments. The mixing induces a coupling of the 
dark Higgs to the Z boson, allowing for a dark Higgs- 
strahlung process e + e~ — > Zh' . Data from the LEP ex- 
periments constrain the parameter £ 2 = (gHZz/ffnzz) 2 ~ 
1CP 2 Ga] for a very light Higgs-like state, which can be 
translated to a bound on the mixing angle 8 < 0.1 which 
is easily satisfied for A^ ~ y\>. Note that this bound is 
conservative and could in principle be weakened consider- 
ably in specific models depending on the decay modes of 
h! . For very light states in the GeV range, Higgs por- 
tal couplings can be probed by rare heavy-flavor me- 
son decays [|46|, potentially down to mixing angles of 
1CP 2 — 10~ 3 , though again, these constraints depend on 
the decays of the GeV-scale dark states. 

Relic Abundance. Assuming that annihilation pro- 
cesses dominate at freezeout, the relic abundance of 
a stable heavy particle is given by J7cm^ 2 — 0.1 X 
((<7u)/pb) _1 , where (av) is the DM annihilation cross 
section. The predicted relic abundance is to be com- 
pared with the observed value from WMAP, Q,DKih 2 = 
0.1123 ± 0.0035 ^. We may classify the annihilation 
channels according to whether the final states belong to 
the dark sector or the SM: 

I: X + X-> VV,Vh',h'h',... 

II : X + X -> V*, h'*, h* -+ i> S M + ^sm 

Case I corresponds to the 'secluded' DM regime p9| , 
in which the relic abundance does not depend on the 
strength of mediation to the SM. The secluded regime is 
typically most important when the DM candidate(s) is 
not the lightest state in the dark sector and direct an- 
nihilation channels are kinematically open in the dark 



sector. For example, if my < m x the annihilation 
cross section is parametrically (av) xx ^vv ~ iro? D /m xl 
so that the observed WMAP abundance requires a D > 
5x 10- 3 (m x /200GeV) g|. 

If no annihilation channels are available directly in 
the dark sector, then it is still possible to obtain the 
proper relic abundance provided the portal operators 
are active and have a sizable coupling. This is case 
II above in which DM may annihilate through the ki- 
netic mixing or Higgs portal operators into the light 
SM final states. As an example, consider the Z^ scalar 
model of Eq. (|). If the DM candidate S is the light- 
est state in the dark sector, it must annihilate into SM 
states if it is to be thermally depleted. One option is 
to utilize the 'direct' Higgs portal operator in Eq. jlq), 
(H^H)(x'X) ^ hSS. In the regime ms "C m/, and away 
from the s-channel resonance at mg = m^/2, the domi- 
nant annihilation mode is into bb. The annihilation cross 
section is (&v) ss _ > . b i ~ iXtm^/irm^ which dictates that 
A x > 0.2 x (m /l /120GeV) 2 . This is a sizable coupling, 
much larger than the bottom quark Yukawa, and so this 
model makes the interesting prediction that the SM Higgs 
decays to DM. In fact in the limit that all other states in 
the dark sector are heavy, this model matches on to the 
minimal Z 2 scalar singlet model of Refs. |t|, [3(j[ Ejj| . 

Next, we consider the relic abundance in multi- 
component models. Typically we might expect that if 
two DM particles have distinct annihilation channels gov- 
erned by different couplings and mass scales, then bar- 
ring some accidental coincidence in these couplings, the 
two DM states will end up with dramatically different 
relic densities and one state will dominantly make up 
the observed cosmological DM. The situation is some- 
what more constrained if the two states have the same 
annihilation modes and symmetry guarantees the equal- 
ity of couplings. This occurs, for example, when two 
heavy states annihilate into light vectors. Given our 
discussion above regarding this channel, we expect then 
if two states Xi > with i = 1,2, have charges Qi and 
masses mi, their relative abundance will be given by 
Oi/r2 2 = {m\/m\)(QyQ\). In this case, a similar relic 
abundance requires a coincidence of masses and charges. 
However, even if a particular species is subdominant in 
the cosmic DM, it still may give signatures in direct and 
indirect detection experiments. We will illustrate this 
below for multi-component direct detection. 

The final interesting possibility we wish to mention 
is the process of semi-annihilation, examined recently 
in Ref. |OJ. Typical annihilation processes change the 
number of DM particles by two units. In contrast semi- 
annihilation reactions change the number by only one 
unit. The simplest example of this occurs in the Z% scalar 
model in Eq. (M) [O . If the dark Higgs h! is lighter than 
the DM xi, then the process XX ~^ X^' wm occur with 
an annihilation cross section (av) xx ^ x h> = 3A 2 /1287rm 2 , 
implying that A > 0.1 x (m x /200GeV) if other annihila- 



tion processes are small. In fact semi-annihilation occurs 
quite generally in multi-component DM models, where 
processes such as XiXj ~^ Xkh' , etc, can take place, and 
we refer the reader to Ref. |ll| for more details. 

Direct Detection. We now discuss the possibility of 
direct detection of DM. Scattering with nuclei is medi- 
ated via a i-channel exchange of vectors or Higgs bosons 
through the portal interactions in (|l|). Vector exchange 
leads to the following DM-nucleon scattering cross sec- 
tion J2l: 
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The most recent limits from the CDMS-II Q and 
XENON100 [|3| experiments have reached a sensitivity 
°f °~ n 5s (3 — 4) x 10~ 44 cm 2 for DM with masses larger 
than around 50 GeV J34J . This translates to a bound 
^9dQ 5s 10 -5 x (rny/lOGeV) 4 for the vector portal. 

Exchange of the SM Higgs boson due the portal inter- 
action also mediates elastic scattering with nuclei. As the 
simplest example consider the Zi scalar model in Eq. (g) . 
The S'-nucleon cross section is (see e.g. [pl| ) 
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where / = ^(n\m q qq\n) /m n ~ 0.35 for where n is 
a proton or neutron |p5| . Constraints from CDMS- 
II H and XENON100 [|| imply A x < 0.03 x 
(m x /50GeV)(m ft /120GeV) 2 . Note that for smaller 
masses in the 10 GeV range, the sensitivity of these ex- 
periments falls off rapidly and larger values of the por- 
tal couplings can be accommodated. Indeed, there has 
been some excitement for DM with masses in the 10 
GeV range, since such a light DM candidate has the 
potential to explain the longstanding DAMA/NAI and 
DAMA/LIBRA pq] signal as well as the recent events 
from COGENT pTJj. This framework can accommodate 
such a light DM candidate as well, but we will not pursue 
this possibility further here. 

Another possibility to explain the DAMA results is in- 
elastic DM H . For a recent study of the status of inelastic 
DM, see Ref. J38). As we have seen, mass splittings may 
occur very easily in particular models with discrete gauge 
symmetries. For example, the models with a Zi discrete 
symmetry naturally lead to inelastic scattering due to the 
off-diagonal coupling to the vector V in Eqs. (Spl), and 
in fact, very similar models were originally considered as 
examples in ||. More generally, inelastic couplings oc- 
cur for a multi-field Zn model when there exists matter 
fields Xn/2 with charge N/2, or two matter fields with 
charges Q and Q + N as in Eq. ([l2]). Small mass split- 
tings of O(100 keV) can lead to a very long-lived excited 
state with a significant relic abundance, as pointed out in 
Ref. p9| , and this leads to the possibility of exothermic 
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FIG. 1: Recoil spectrum of two-component scattering at 
XENON100. Here we consider a Z 6 model in Eq. (jl|) 
with two DM particles S D X3 an d X2- We have assumed 
masses for the DM particles ms — 3 TeV with a splitting 
of S = mp — ms — 150 keV leading to inelastic transitions, 
and m X2 — 30 GeV. Other parameters entering into the relic 
abundance and scattering cross section are k = 8 x 10~ 4 , 
qd ~ 0.01 and mv = 10 GeV. This point is consistent with 
current direct detection constraints [B2, B3] and should lead 
to ~ 15 events at XENON100 assuming an effective expo- 
sure of 3000 kg-days. The dashed(dotted) lines show the 
elastic(inelastic) components while the solid line shows the 
total rate for summer (June 2nd). To compute the rate, we 
have used a standard Maxwellian velocity distribution with 
v esc — 500 km/s. 



down-scattering with nuclei. Recently, this has been uti- 
lized in Refs. EQ, [ll[ as a possible explanation of DAMA. 

There is also the possibility of 'semi-elastic' scattering, 
investigated recently in Ref. [j42|| , which has a very dis- 
tinctive recoil spectrum. Semi-elastic scattering can oc- 
cur due to the exchange of a vector with an off-diagonal 
coupling and a Higgs with a diagonal coupling, as in p2| , 
or through the exchange of a vector with both diagonal 
and off-diagonal couplings as in Eq. (|l3|). In the latter 
case, in order to have the small mass splitting required for 
an inelastic transition, the bare masses of the two fields 
must be approximately equal, and we have not provided 
a symmetry reason for this here though in principle this 
could be addressed with further model building. 

Finally, we turn to the direct detection signatures of 
multi-component DM, which has been considered previ- 
ously in Ref. jl3| . Naively, one might expect that detect- 
ing two components is only feasible if both DM particles 
have similar relic densities. As we discussed above, this is 
not generally true unless the both particles have similar 
masses and annihilation channels. However, the rate de- 
pends both on the relic density and the scattering cross 
section. A particle with a large annihilation cross section 
that is efficiently depleted early in the universe might also 
be expected to have a larger scattering cross section with 
nuclei. Furthermore, the particle with the dominant relic 
density can scatter inelastically with nuclei, which leads 



to a smaller cross section compared to the elastic case. 
These two facts suggest that it may indeed be possible 
to have significant rates for both DM particles at direct 
detection experiments. 

As an example, consider the Z§ model in Eq. (IIG) 
with two stable particles, the real component S of X3 
and X2- These states are stable due to the presence of 
Z3 and Zi symmetries as explained above. Assuming 
annihilations to dark gauge bosons V, whichever parti- 
cle is heavier will dominantly comprise the cosmological 
DM, and we assume this to be S. Because the S is split 
from its partner P, an off-diagonal gauge interaction as 
in Eq. (|J) mediates inelastic scattering with nuclei, and 
this will be quite suppressed for splittings of 0(100 keV). 
On the other hand \i wm undergo elastic transitions, 
and because of its dilute presence in the galactic halo, 
we are free to increase k to much larger values without 
conflicting with direct detection constraints. If there is a 
large mass separation between the two DM components, 
a very distinctive recoil spectrum can be observed at cur- 
rent and next generation experiments, as we illustrate in 
Fig. I 

Collider Signatures. When considering new particles 
from a dark sector, because of the weak coupling to the 
SM it is typically hopeless to directly produce such states 
at colliders. However, an interesting opportunity arises if 
a hidden valley scenario is operative [[J7| , in which case we 
first produce new heavy 'connector' states with SM quan- 
tum numbers in abundance, which subsequently decay 
into the dark sector states. While there are many pos- 
sibilities for connector states around or above the weak 
scale, we will focus on the minimal case of the SM Higgs 
boson as our connector. The SM Higgs is in fact an ideal 
choice for a bridge to the dark sector because of its sizable 
production cross section, cr^^° ~ 50(10) pb for LHC at 
y/a = 14(7) TeV, and narrow width. For m h < 130 GeV, 
the Higgs decay width is governed by the small bottom 
quark Yukawa yb ~ 1/60, so that a light SM Higgs is 
quite susceptible to new decay modes Q. We will uti- 
lize this property of the Higgs to gain access to states in 
the dark sector. 

From Eq. ( |l8| ) we see there are two possible operators 
that will allow for new decays of the SM Higgs. The 
operator <fvcf)WH allows for the decay of the SM Higgs 
to dark Higgs and gauge bosons, with partial widths 
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where x, = rrii/mh. The SM Higgs couples to the 
dark vectors due to the mass mixing induced by C por t a l ■ 
Notice that in the limit Xh',v "^ 1) the branching to 
dark Higgs bosons and dark vectors is equal, T h ^ h > h > — 
T/i^yy — A 2 i w 2 /87TTO/ l , a simple consequence of the 
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FIG. 2: SM Higgs branching ratios. Here we are 
considering the Z2 scalar model of Eq. (pi), with a 
spectrum (mv,ins,iTip,mh') — (20, 30, 50, 60) GeV, and 
their corresponding branching ratios are indicated by 
(VV,SS,PP,h'h') = (blue, red, orange, green). We have 
taken portal couplings to be A0 ~ yb — 1/60, and A x = 1/20. 
The SM 66, ZZ and WW modes are indicated by the (solid, 
dotted, dashed) black lines. Below the WW threshold, the 
SM Higgs decays dominantly to hidden sector states. 



Goldstone-equivalence theorem since the Higgs portal 
Jig) leads to identical couplings of the SM Higgs to the 
radial and Goldstone modes of the dark Higgs. We should 
compare these partial widths with the partial width to 
a pair of bottom quarks, T h _^. b i ~ 3y 2 m/i/167T, telling 
us that for A^ ~ yb the branching into dark Higgs and 
gauge bosons will be significant. Note that this value of 
A^ corresponds to a mixing angle 6h ~ 10 -3 — 10~ 2 , well 
below any existing constraints. 

If the theory contains scalar matter fields %, where \ 
now generically denotes either the stable DM state or 
an unstable Zm partner, then the second Higgs portal 
operator x^xH^H in Eq. (|18| ) allows a decay h 
with partial width 
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A similar conclusion applies: as long as A x is comparable 
to yb, the SM Higgs will have a sizable branching into 
DM or its Zn partners. We show the branching ratios of 
the SM Higgs in Fig. Bfor a sample model and spectrum. 

After the primary decay of the SM Higgs into either 
h'h', VV, or x^X> the observed signature depends on the 
cascade decays of the dark states. This becomes highly 
dependent on both the field content and the spectrum 
in the dark sector. Even a small number of additional 
unstable Zn matter states in the dark sector can lead to 
very complicated signatures due to long cascades, partic- 
ularly when there is mass mixing or nontrivial Zjq inter- 
actions as in (|l4|) . It is therefore useful to first consider a 



simplified situation where there is one stable DM x an d 
a heavier unstable Zjq partner x* i in which case we can 
classify the possible decay patterns as follows: 

• x stable (^), 

• X* -> 2 X> V\, hx , 

. v^2 X ,2 X *,xx*,i + e-,jj , 

• h'^2 X ,2 X *,xX*,VV,jj . 

Note that this is schematic, since for the matter states 
X, X* we have not distinguished between particle and an- 
tiparticle, and furthermore in particular models the real 
and imaginary components of matter fields may receive 
mass splittings. We observe two important points from 
the decay patterns above: 1) It is quite common to de- 
cay into a stable DM state, which will lead to a missing 
energy signature, and 2) dark states may cascade decay 
back to the SM, leading to a dilepton or dijet signature. 
Thus the generic signatures in this class of models are 
multi-lepton, multi-jet final states with missing energy. 

While the origin of the missing energy signature is 
fairly obvious, it is useful to examine in more detail how 
it is possible for dark sector states to decay back to the 
SM. This occurs due to the kinetic mixing and Higgs 
portals in Eq. ([lq). For the vector, the kinetic mixing 
operator V ilv B^ v leads to the coupling fli~9| ) of V to the 
SM fermions. The partial width for a vector to decay to 
a generic SM fermion / is T v ^fj ~ (N c /3)an 2 my. If 
there are no direct channels open to dark sector states, 
the vector will decay to SM fermions. For very light 
vectors, below the hadronic threshold, the vectors de- 
cay into e + e~ pairs or, if kinematically allowed, /x + // _ 
pairs. Above the dipion threshold, the vector can decay 
to hadrons; the hadronic yield is governed by the well- 
measured form factor R = (J e+e -^ hadron Jo- e + e -^^+ M - 
for vectors lighter than several GeV ^4], H. At higher 
masses, one can use the quark description to calculate 
the decay to hadrons |30| . In summary, the vector decays 
fairly democratically, and in particular there is always a 
significant 0(1) branching to dilepton final states. Fi- 
nally, the decays are typically prompt for n > 10 -5 , but 
there may be displaced vertices for smaller kinetic mixing 
parameters. Similarly, the dark Higgs h! can decay via 
the Higgs portal to the heaviest kinematically available 
SM fermion pair. Typically this is the 6-quark pair, and 
so the partial decay is given by T h ,^ b i ~ 30| 1 ?/ 2 m/ v /167r, 
which is a prompt decay for Oh ^ 10 -4 . A leptonic sig- 
nature from the Higgs portal is not possible unless the 
dark Higgs is very light, in the hundreds of MeV range, 
in which case there is some tension with technical natu- 
ralness. From now on we will focus on the clean leptonic 
signature via the kinetic mixing portal. The possibility 
of the SM Higgs decaying to multi-lepton final states in 
this manner was discussed in Ref. ptj . 



Because of the kinetic mixing suppression in the decay 
V —> ff, if kinematically allowed the dark vector will 
prefer to decay directly into states in the dark sector. 
However, this does not necessarily imply an absence of 
SM particles in the final states, since other dark sector 
states may decay back to the SM through an off-shell V* . 
This occurs when mass splittings or mixings lead to off- 
diagonal gauge interactions, as given e.g. in Eqs. (^,p|,^3[). 
As an illustration, consider the Z2 scalar model with the 
gauge interaction in Eq. (0), and a spectrum such that 
my 3> mp > ms, so that on-shell vectors decay domi- 
nantly via V — > PS. The heavy scalar partner will decay 
via the three-body process P — > SV* — > Sff, and in the 
limit of large my and small fermions mass the partial 
width is given by 
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where /(x) = 1 — 8x — 12a; 2 log: 
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The quantity 



(g v + g\) can be extracted from Eq. ( |19| ) and depends on 
the ratio my/mz', light vectors my/mz <C 1 couple to 
electric charge, (g v + g A ) — C%Q% M , while heavy vectors 
my/mz >• 1 couple to hypercharge, (g v + g A ) ~ (Y 2 + 
Yft)/2. For a particular fermion final state to be allowed 
kinematically, the mass splitting must be large enough, 
mp — ms > 2m f. The lifetime depends sensitively on 
the masses of the vector and Z^ partner P as well as the 
kinetic mixing, but for masses in the 10-100 GeV range 
the displaced vertices may occur even for large values of 
the kinetic mxing n < 10 -3 . 

An interesting probe of multi-component Zjy models 
would be the presence of two or more distinct stable 
particles in the cascade decays of the SM Higgs. This 
may happen from an initial asymmetric decay of the SM 
Higgs into two distinct states, h —> XiXj with i ^ j which 
subsequently cascade into distinct stable states plus SM 
fermions. Even if the initial decay is symmetric, e.g. 
h —> VV, if V has comparable branchings into distinct 
final states then there will be events which contain differ- 
ent stable particles. Consider as an example the Z4 model 
of Eq. (||) with two stable candidates S and xi- Provided 
the vector is much heavier than the DM candidates, the 
vector will decay via V — >• SP and V — > XiXi with partial 
widths proportional to the square of the charges. Thus, 
the yield of SP vs. X1X1 w iH be 4:1. The signature in 
this case would be £ + l~ f} with missing energy coming 
from two S particles on one side of the event and two Xi 
particles on the other. There are of course many other 
possibilities and in general it would be worthwhile to in- 
vestigate the extent to which the properties (mass, spin, 
etc.) of the different stable states can be determined. 
There has been recent work in this direction in Ref. pM . 
Along these lines, an interesting question is whether col- 
liders have something to say regarding the nature of the 
discrete symmetry at hand, i.e. can one distinguish be- 
tween, say, Z2 and Z3 symmetries? This question has 



received attention recently in Ref. [ p2| , where the au- 
thors study the difference in the possible decays of parent 
particles obeying a Z2 vs. Z3 symmetry and how such 
differences manifest in the invariant mass distributions 
of the outgoing SM particles. Another avenue of deter- 
mining the underlying discrete gauge symmetry in these 
models is to extract the charges of the Higgs boson and 
matter fields from collider data. Measuring specific de- 
cay widths and branching ratios of e.g. the gauge boson 
into matter states can provide some information in this 
regard. 

If the dark sector is very light compared to m^, say 
at a scale of O(GeV) or less, the decay products will 
be highly boosted, leading to a novel "lepton jet" signa- 
ture [ p3| : multiple highly-collimated groups of leptons in 
the final state. In fact, it has recently been suggested 
that a SM Higgs decaying to lepton jets could have been 
missed at the LEP experiments [Q, allowing for a re- 
laxed bound on the SM Higgs mass. As in the unboosted 
case, various aspects of the collider phenomenology, such 
as event topology, lepton multiplicity, and missing en- 
ergy, is highly dependent on the spectrum in the dark 
sector. Many of the suggested search strategies in [ pi] , 
and past searches at LEP and Tevatron will be relevant 
for the class of models considered here. In particular, 
if we restrict to a Higgs mass above the LEP bound, 
nih > 114 GeV [J45J, searches at the Tevatron for 'dark 
photons' H), NMSSM pseudoscalars H, trileptons J57), 



and same-sign dileptons |p8| have the potential to con- 
strain particular models and spectra. 

To summarize, if moderate couplings to the portals 
on the order of the bottom quark Yukawa exist, there is 
an opportunity to probe the dark sector at high energy 
colliders via the decays of the SM Higgs boson. Look- 
ing forward, it seems useful to proceed in two directions 
due to the model dependence of the signatures. First, 
it would be worthwhile to study particular benchmark 
models containing some of the more complex signatures 
of multi- lepton, multi-jet final states with missing energy. 
Complementary to this, the parameter space of the more 
minimal models, particularly those with a single matter 
field, can be thoroughly explored. 

Discussion. Gauge symmetry provides a plausible 
origin for the discrete symmetries that can stabilize DM. 
In this paper we have surveyed a class of DM mod- 
els based on an abelian Zjy discrete gauge symmetry. 
We have investigated models with minimal field content, 
multi-field Zjq models, and models with chiral matter. 
We have found that Zn symmetries may lead to multi- 
ple stable DM candidates when N is not prime. The dark 
sector may couple to the SM through the kinetic mixing 
and Higgs portals, and we have given an outline of the 
basic phenomenology of this scenario, highlighting some 
novel direct detection and collider signatures. In partic- 
ular, the SM Higgs boson may provide a portal through 
its decays to new states in the dark sector. 



Our focus in this work has been on the gauge origin of 
the discrete symmetries and its potential impact on DM 
physics. There are a number of avenues for future work. 
Regarding indirect probes of these models, DM annihila- 
tion in the galaxy and in the sun can lead to novel astro- 
physical signatures. As is well-known by now, a very light 
gauge boson provides a 'dark force' [M that can poten- 
tially explain the cosmic ray signatures of PAMELA |22 ] 
and FGST p3[. In this context, there is scope for the 
exploration of the low energy phenomenology of these 
models in fixed target p6|, |27| and e + e~ experiments [ p4[ . 
Also, DM annihilation in the sun to light metastable me- 
diators, such as dark Higgs or gauge bosons, can yield 
electromagnetic signatures in gamma-ray telescopes |59(] 
such as FGST [23]. In terms of cosmology, as recently 
pointed out in [3(J GeV-scale states from a dark sector 
can help to resolve the lithium problem of Big Bang Nu- 
cleosynthesis. It would be valuable to study these and 
other connections with cosmology and astrophysics for 
the DM models presented here. 

In terms of model-building, we have not attempted to 
address the electroweak hierarchy problem, or the addi- 
tional naturalness issues that come with scalars in the 
dark sector. In this regard, it would be interesting to 
adapt these models to a supersymmetric framework along 
the lines of |53|, |6l| or a warped-extra-dimensional frame- 
work []62[ . This can dramatically alter the collider phe- 
nomenology, as now one may use the supersymmetric or 
Kaluza-Klein states as connectors to the dark sector. 

We have focused on DM in a hidden sector, but one 
could also consider new U(l) gauge symmetries under 
which the SM is charged. There are a few anomaly-free 
global U(l) symmetries that one may gauge within the 
SM: B—L (with right-handed neutrinos) or a difference in 
lepton numbers, e.g. L e — L^ J63j. However, more general 
U(l) symmetries can be gauged provided new fermion 
representations are added to cancel anomalies, and these 
could be used to obtain discrete gauge symmetries. 

Stability (at least on cosmological time scales) is one of 
the few known properties of DM and hints at a new sym- 
metry in the dark sector. Should experimental evidence 
for DM be found, new stabilization symmetries may be 
required to understand the data. In the meantime, ex- 
ploring new possibilities for these symmetries may sug- 
gest new phenomena associated with DM and therefore 
new search strategies. Our investigation here illustrates 
this basic observation, as we have been led by symmetry 
considerations to predict new states, new interactions, 
and ultimately new signatures for direct detection and 
collider experiments. 
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